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,— 1 Abstract 

in 

We use perturbation theory to construct perfect lattice actions for 
quarks and gluons. The renormalized trajectory for free massive quarks 
is identified by blocking directly from the continuum. We tune a pa- 
1} rameter in the renormalization group transformation such that for 1-d 

configurations the perfect action reduces to the nearest neighbor Wilson 
fermion action. The fixed point action for free gluons is also obtained by 
blocking from the continuum. For 2-d configurations it reduces to the 
standard plaquette action. Classically perfect quark and gluon fields, 
quark-gluon composite operators and vector and axial vector currents 
are constructed as well. Also the quark- antiquark potential is derived 
from the classically perfect Polyakov loop. The quark-gluon and 3-gluon 
perfect vertex functions are determined to leading order in the gauge 
coupling. We also construct a new block factor n renormalization group 
transformation for QCD that allows to extend our results beyond per- 
turbation theory. For weak fields it leads to the same perfect action as 
blocking from the continuum. For arbitrarily strong 2-d Abelian gauge 
fields the Manton plaquette action is classically perfect for this transfor- 
mation. 
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1 Introduction 



Cut-off effects are the major source of systematic errors in numerical simulations 
of lattice QCD. For the standard Wilson action these effects are of the order of 
the lattice spacing. Hence they vanish rather slowly as the continuum limit is ap- 
proached. Symanzik's perturbative improvement program systematically eliminates 
cut-off effects by introducing irrelevant higher- dimensional operators in the lattice 
action [1]. For QCD this program was realized by Sheikholeslami and Wohlert [2]. 
Their so-called 'clover action' is now often used in numerical simulations. A related 
attempt to damp cut-off effects uses the concept of tadpole improvement [3]. 

Recently Hasenfratz and Niedermayer proposed a new approach to eliminate cut- 
off effects by using nonperturbative approximations of 'perfect actions' [4]. By defi- 
nition a perfect action is free of any cut-off artifacts. It is located on a renormalized 
trajectory emanating from a fixed point on the critical surface. The crucial obser- 
vation of [4] is that the construction of the fixed point action in an asymptotically 
free theory is a classical field theory problem that can be solved nonperturbatively. 
In fact, the fixed point action is a classically perfect action. Using this action even 
in the quantum theory Hasenfratz and Niedermayer observed that cut-off effects are 
practically eliminated in the 2-d 0(3) nonlinear cr-model at correlation lengths as 
small as 3 lattice spacings. The action has been combined with the classically perfect 
topological charge in a study of the scaling properties of the topological susceptibil- 
ity [5]. Recently the renormalized trajectory for a ^-function renormalization group 
transformation was investigated for the 0(N) model in [6]. Wieczerkowski and Xy- 
lander combined the concepts of perfect and Symanzik improved lattice actions in 
a study of the 2-d hierarchical 0(N) nonlinear cr-model [7]. For the Gross-Neveu 
model in the large N limit it has been shown analytically that the fixed point ac- 
tion is exactly perfect at the quantum level even at arbitrarily small correlation 
length [8] . Also it has been realized that a classically perfect action is automatically 
1-loop quantum perfect [9]. Finally the perfect action approach has been applied 
successfully to pure SU(3) gauge theory [10]. Other constructions of perfect lattice 
actions for gauge theories are contained in [11]. A more formal discussion of block 
transformations for gauge fields is given in [12], and early numerical work on SU(2) 
can be found in [13]. 

In the present paper we concentrate on the inclusion of the quarks. Using per- 
turbation theory we construct a lattice QCD action that is perfect to first order in 
the gauge coupling for arbitrary fermion mass. This is an essential step towards a 
fully nonperturbative perfect action for QCD. Since nonperturbative effects are not 
yet included our action is not intended to be used in simulations of light quarks. For 
heavy quarks, however, the leading cut-off effects are related to the heavy fermion 
mass. In our approach these effects are completely eliminated to first order in the 
gauge coupling. Since the gauge coupling is weak at the scale of a heavy fermion 
mass, our perturbatively perfect action is promising for numerical simulations of 
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heavy quarks. The inclusion of nonperturbative effects will be discussed in a sep- 
arate publication [14]. Then we will be able to treat light quarks as well. The 
question of locality of the action plays an important role throughout this paper. We 
distinguish three types of locality. An action is called 'ultralocal' if its couplings are 
limited to a finite set of neighbors. It is called 'local' if the couplings decay expo- 
nentially at large distances, and 'nonlocal' if they decay only powerlike. The main 
technical tool of this paper is a renormalization group transformation that defines 
lattice fields by blocking directly from the continuum. In section 2 we discuss this 
transformation for free massive quarks. The perfect action is determined and a pa- 
rameter in the renormalization group transformation is tuned such that the action 
becomes ultralocal for 1-d configurations. In section 3 we apply the technique of 
blocking from the continuum to free gluons. Again the parameters in the renormal- 
ization group transformation are optimized for locality of the fixed point action. For 
2-d configurations the action reduces to the standard ultralocal plaquette action. In 
section 4 classically perfect quark and gluon fields as well as composite operators 
and currents are constructed. Also the static quark- antiquark potential is derived 
from the classically perfect Polyakov loop, and it is compared with the potentials 
of the standard Polyakov loop and of the continuum theory. In section 5 quarks 
and gluons are coupled, and perfect quark-gluon and 3-gluon vertex functions are 
derived to first order in the gauge coupling. Section 6 contains the construction of a 
nonperturbative block factor n renormalization group transformation that is (upon 
iteration) equivalent to blocking from the continuum for weak fields. This enables 
the extension of our results beyond perturbation theory. The new transformation 
has attractive features. In particular, for arbitrarily strong 2-d Abelian gauge fields 
the Manton plaquette action [15] is classically perfect. Finally section 7 contains 
our conclusions. A synopsis of this paper has been written up in [16]. 



2 Blocking free quarks from the continuum 

Block variable renormalization group transformations with a finite blocking factor 
are well established in lattice field theory and statistical mechanics. Iterating such 
transformations leads — for a suitable choice of parameters — to fixed points of 
the renormalization group, which are directly related to the continuum limit [17]. 
It has been discovered a long time ago (but has not been very well-known) that 
one can reach these fixed points also in one step by blocking directly from the 
continuum. Wilson has applied this method to scalar fields and mentioned possible 
generalizations [18]. A formal discussion is given in [19]. Cronjager and Mai applied 
it to free fermions [20]. It has also been used in the nonlinear cr-model [4, 21]. Finally, 
in a recent perturbative construction of lattice chiral fermions we have derived the 
perfect action of fermions and gauge fields by blocking from the continuum [22] . 

Let us consider a quark field ijj^ijj in the continuum. We construct a lattice quark 
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field by averaging the continuum field over hypercubes c x centered at the points x 
of a d- dimensional Euclidean lattice 

V x = f d d yi;(y), V x = j d d y ^{y). (2.1) 

Jc x J c x 

As a rule, in this paper we use capital letters for lattice fields and small letters for 
continuum fields. The hypercubes have unit length and hence the lattice spacing 
is 1. All physical quantities in this paper — also those of the continuum theory — 
are measured in lattice units. Integrating out the continuum degrees of freedom one 
arrives at an effective action for the lattice fields 

exp(-S[tf,tf]) = / ^exp(- S [^^])[I^(i- / d?yj,{y))8{* x - f d d y^{y)). 

(2-2) 

Here s\^)^\ is the continuum action for free quarks. In momentum space it takes 
the form 

*hM] = j^y J d d p ^(-p)^(p)"V(p), S f (p) = (i^p.+m)- 1 , (2.3) 

where 8^{jp) is the free quark propagator in the continuum. The lattice fermion field 
is given by 

d a 

*(p) = Yl + 27r/ )n(p + 2 *0, n (p) = II -> = 2 sin (^/2)- (2-4) 

l£Z d f =1 Vil 

Note that the momentum of the lattice field is restricted to the Brillouin zone B = 
] — 7r, iv] d . Due to the summation over the integer vector I the lattice field is periodic 
over B. To perform the integration over the continuum fields we go to momentum 
space and represent the Grassmann ^-function of eq.(2.2) as an integral over auxiliary 
lattice fermion fields fj, rj 

exp(-5[*,*]) = J V^VfjVrjexp { - J d d p ^(-p)^(p)"V(p)} 

x exp { — !— / d d p \m-p) - V i>{-p - 27r/)n(p + 2irl)] V (p) 

7/(-p)[*(p) - V>(? + 27r/)n(p + 2tt/)] + arj(-p) V (p)\ }. 



lez" 

; ; ; 

(2.5) 



The ^-function renormalization group transformation of eq.(2.2) corresponds to a = 
0. This transformation is chirally invariant and its fixed point (at m = 0) describes 
free lattice chiral fermions [23]. The corresponding fixed point action is nonlocal. 
Cronjager has shown that its couplings decay as r 1 ~ d at large distances r [20]. 
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Hence it is consistent with the Nielsen- Ninomiya no-go theorem. In [22] we have 
shown that the nonlocality does not cause problems even in the interacting chiral 
gauge theory. In particular, the correct axial anomaly is reproduced on the lattice. 
Here we are interested in QCD which is a vector-like theory. Therefore we allow 
chiral breaking terms in the renormalization group transformation and hence in the 
fixed point action. In principle one may allow such terms also in a chiral gauge 
theory. 1 This possibility was explored in [16]. In fact, for numerical investigations 
it is essential to have a local action. Therefore in eq.(2.5) we have introduced a 
chiral symmetry breaking mass term arjr] for the auxiliary field which smears the 
^-functions of eq.(2.2) to a Gaussian distribution of width a. For a > the perfect 
action turns out to be local [23] . Later we will see explicitly that the chiral breaking 
introduced in the renormalization group step does not affect the spectrum of the 
theory. 

Next we combine the integral over the Brillouin zone in eq.(2.5) with the sum over 
integer vectors I to an integral over the entire momentum space of the continuum 
theory 



Performing the Gaussian integral is equivalent to solving a classical equation of 
motion for the continuum quark field. The solution of this equation is 



iez d 

Performing the Gaussian integral over the auxiliary field is again equivalent to solv- 
ing a classical equation of motion. In this case the solution is 




(2.6) 




(2.7) 




(2.8) 



rj c (-p) = -y(- P )A f ( P )-\ Vc (p) = -A'(p) _1 *(p) 



(2.9) 



We thank P. Hasenfratz for pointing this out to us. 
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with 

A f (p) = 6 f {p + 27r/)n(p + 2tt/) 2 + a, (2.10) 

iez d 

such that together with eq.(2.7) we obtain 

M-p) = ^(-pJA^p)-^^)^^), MP) = 6 f (p)IL(p)Af(p)- 1 ^(p). (2.11) 

These relations between the continuum and lattice quark fields will be useful later 
when we construct composite operators and determine the perfect quark-gluon ver- 
tex function. Inserting eq.(2.9) back into the exponent of eq.(2.8) one obtains the 
perfect action 

= -±- d f d d P ^(-pjA'Cp)- 1 *^), 

l Z7r J Jb 

a/ (p) = ^[^(^ + 27r/ ") + m ]" llI ( ?5 + 27r/ ) 2 + a - ( 2 - 12 ) 

iez d 

The sum in the lattice fermion propagator A-^(p) converges in any dimension because 
II(p + 27r/) 2 suppresses contributions from large [l^. The lattice action describes a 
complete renormalized trajectory parametrized by the mass m. At m = the 
trajectory intersects the critical surface in the fixed point. 

The introduction of a > leads to a local perfect action with explicitly broken 
chiral symmetry even at m = 0. We now vary a such that the action becomes 
as local as possible. It turns out to be sufficient to optimize the locality for 1-d 
configurations. In one dimension the sum in eq.(2.12) can be performed and the 
fermion propagator takes the form 

A / (p) = — - icoth(m/2) - i cot(p/2)]" 1 + a. (2.13) 

Choosing a = (m — m)/m 2 , where m = e m — 1, the propagator reduces to 

A f (p)= (—Yiismp + m + lf)- 1 . (2.14) 

This corresponds to the standard nearest neighbor Wilson fermion action. Hence the 
perfect action is ultralocal for 1-dimensional configurations. At m = the optimal 
parameter is a = 1/2. The same fixed point was identified before by numerically 
optimizing the locality in two dimensions [23]. It turns out that the perfect action 
is still very local even in d = 4. To see this we go back to coordinate space and we 
write the action as 

s[*, *] = X) **[7^(* -y) + A(» - y)]*v (2-15) 



x,y 



Performing a Fourier transform of the inverse fermion propagator for the optimal 
value a = (m — m)/m 2 we find p^z) and A(z). Their values for m = 0, 1 and 2 are 
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Z\ 


^2 


Z 3 


z 4 


m = 


m = 1 


m = 2 


1 











0.1400 


0.05503 


0.01856 


1 


1 








0.0314 


0.01088 


0.00314 


1 


1 


1 





0.0096 


0.00302 


0.00076 


1 


1 


1 


1 


0.0033 


0.00100 


0.00023 


1 


2 








0.0015 


0.00024 


0.00003 


1 


2 


1 





0.0006 


0.00010 


0.00001 


1 


2 


1 


1 


0.0002 


0.00003 


0.00000 


1 


3 








0.0001 


0.00001 


0.00000 


2 











0.0040 


0.00054 


0.00003 


2 


1 








0.0012 


0.00018 


0.00002 


2 


1 


1 





0.0002 


0.00004 


0.00001 


2 


1 


1 


1 


-0.0002 


-0.00002 


0.00000 


2 


2 








-0.0001 


-0.00003 


-0.00001 


2 


2 


1 





-0.0001 


-0.00003 


0.00000 


2 


2 


1 


1 


-0.0002 


-0.00002 


0.00000 


3 











0.0004 


0.00004 


0.00000 


3 


1 








0.0001 


0.00001 


0.00000 



Table 1: Some values of p\(z) in d = 4 for masses m = 0, 1 and 2. The parameter 
a is chosen such that the action is ultralocal for 1-d configurations. Note that p\(z) 
is symmetric under permutations of z 2 , z 3 and z 4 . 



7 





%2 


%3 


z 4 


m = 


m = 1 


m = 2 














1.8530 


1.26893 


0.84425 


1 











-0.0616 


-0.03039 


-0.01206 


1 


1 








-0.0290 


-0.01067 


-0.00325 


1 


1 


1 





-0.0146 


-0.00446 


-0.00111 


1 


1 


1 


1 


-0.0077 


-0.00205 


-0.00043 


2 











0.0010 


0.00034 


0.00009 


2 


1 








-0.0007 


-0.00012 


-0.00001 


2 


1 


1 





-0.0007 


-0.00013 


-0.00002 


2 


1 


1 


1 


-0.0005 


-0.00008 


-0.00001 


2 


2 








0.0004 


0.00008 


0.00001 


2 


2 


1 





0.0002 


0.00004 


0.00001 


2 


2 


1 


1 


0.0001 


0.00002 


0.00000 


2 


2 


2 





0.0001 


0.00002 


0.00000 


2 


2 


2 


1 


0.0001 


0.00001 


0.00000 



Table 2: Some values of X(z) in d = 4 for masses m = 0, 1 and 2. Again the 
parameter a is chosen such that the action is ultralocal for 1-d configurations. Note 
that X(z) is symmetric under permutations of the z^. 

given in table 1 and 2, and are illustrated in fig.l and 2. For a lattice spacing of 0.3 
fm (which may be realistic for a simulation with a perfect action) the value m = 2 
roughly corresponds to the mass of the charm quark. 

Let us discuss the spectrum of the fermionic lattice theory. We investigate the 
2-point correlation function of a fermionic operator with definite spatial momentum 

pe]-7r,7r] rf -\ 

i r 

*(pk = Tr/ d Pd y(p)exp(i Pd x d ), (2.16) 

that creates a fermion at Euclidean time x d and the operator Vf^o) that annihilates 
a fermion at the point x at Euclidean time 0. The correlation function is given by 

(^(^o)*^^} = -!- / dp d TTA f (p)exp(ip d x d ) 

2V J -IT 

= — I dp d Tr [ V S f (p + 27r/)n(p + 2tt/) 2 + a] exp(ip d x d ) 
= — / dv d y 



, (p + 27T/) 2 +p2 +TO 2 
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Figure 1: Exponential decay of \p ll {i,i,i,i)\ (along the diagonal) for masses m = 
(diamonds), m = 1 (crosses) and m = 2 (squares). The lines are least square fits of 
\p ll (i,i,i,i)\ oc exp(— c(m)i) with c(0) = 4.939, c(l) = 5.857 and c(2) = 6.969. 

= C(p + 2Tvl)exp(-E(p + 2Tvl)x d ) + a6 Xdfi . (2.17) 

fez**- 1 

The sum over Id has been combined with the integral of p& over ] — 7r,7r] to an 
integral over R, characteristic for the continuum theory. The sum over spatial / 
leads to infinitely many poles of the integrand, and thus to infinitely many states 
that contribute an exponential to the 2-point function. The energies of these states 
are given by the location of the poles, E(p + 2ttI) = —ipd, with 

E(p + 2tvI) 2 = -p\ = (p + 2tvI) 2 + m 2 . (2.18) 

This is exactly the energy of a particle with mass m and momentum p + 2ttI, i.e. 
the spectrum of the lattice theory is identical with the one of the continuum — 
cut-off effects are completely eliminated. This is possible only because the lattice 
momentum p, which is restricted to the Brillouin zone, combines naturally with 
the integer vectors / to cover the whole momentum space of the continuum theory. 
Therefore the continuous Poincare symmetry is restored in the spectrum, although 
it is explicitly broken in the action. 

The nonperturbative construction of an approximately perfect action for QCD 
will require some numerical work. Then it is essential to use a renormalization group 
transformation with a finite blocking factor. In section 6 we will define a factor n 
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Figure 2: Exponential decay of (along the diagonal) for masses m = 

(diamonds), m = 1 (crosses) andm = 2 (squares). The lines are fits of \X(i,i,i,i)\ oc 
exp(— c(m)i) excluding the point (0,0,0,0). Here c(0) = 4.967, c(l) = 5.850 and 
c(2) = 6.942 are very close to the corresponding values for \p ll {i,i,i,i)\. 



transformation whose perfect action is identical with the one of blocking from the 
continuum at least for weak fields. To optimize the parameters of the factor n 
transformation we want to relate them to those of blocking from the continuum. In 
the discrete transformation the fermionic variables ty x , ty x of an n d block on the fine 
lattice are averaged to form coarse lattice variables 

*i' = ^£*-' *i' = ^£*- (2-19) 
n a ^— ' n a ^— ' 

x£x' x£x' 

located at the block centers x' that form a lattice of spacing n. (By x £ x' we 
denote that the fine lattice point x belongs to the block x' .) The factor b n is a 
wave function renormalization that is necessary because on the coarse lattice we 
again work in lattice units (now of the coarse lattice). The renormalization group 
transformation with block factor n takes the form 

exp(- S' [*',*']) = J V^V^VfjV V exp ^ - j^y J d d p *(-p)A / (p)- 1 *(p)} 

X exp | J ^ d d p \a n T){-p)fi{p) 

+ [*'(-p)-6 B ^*(-p-^)n„(p + — )]^(p) 
* — ' n n 
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+ ^(-p)[*'(p)-6„X;*(p + — )n»(p + — )]1}. (2-20) 

/' 

Here 5' =] — iv /n,iv /n] d is the Brillouin zone of the coarse lattice. As before we 
have introduced a smearing parameter a n . The sum over /' now extends over integer 
vectors with components restricted to V £ {l,2,...,n} and II„(p) = II (rap)/ II (p). 
Integrating out the fine lattice and auxiliary fields we obtain the propagator of the 
coarse lattice field 

n a — ' n n 

v 



4 £ [ £ + M,) + m]-n(£±^ + 2W)» + a 

n a — ' L * — ' n n 

v iez d 

x n„(^ ) 2 + a ri 



n 

l'eZ d I' 
bl , _ _ _ hi 



n d=l X ['^(Pai + 2ttQ + nrnj-^d) + 27r/') 2 + -£a + o„. (2.21) 

This is indeed the perfect propagator for blocking from the continuum provided 
that b n = 7J,( rf_1 )/ 2 and a/n + a n = a. The value of b n is understood on dimensional 
grounds because the dimension of a fermion field is (d — l)/2. The mass is a factor 
n larger because it is now measured in units of the coarse lattice. One obtains [24] 

a n = a(l-~y (2.22) 

This means that e.g. the optimal value a 2 for a factor 2 renormalization group 
transformation is half the value of a that is optimal for blocking from the continuum. 
Also = a because blocking from the continuum is equivalent to blocking with an 
infinite blocking factor. 



3 Blocking free gluons from the continuum 

In this section we apply the technique of blocking from the continuum to free gluons. 
As long as interactions are switched off the gluons behave like decoupled Abelian 
gauge fields. Therefore we suppress color indices and denote the gluon field in the 
continuum by a^. In analogy to the fermions we want to average the continuum 
gauge field to obtain a blocked lattice field. Lattice gauge fields naturally live on 
the links connecting neighboring lattice points. The lattice field is therefore obtained 
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by an integration over the two hypercubes centered at these lattice points. We define 
a noncompact lattice gauge field 

A-»,x = / d d y (l + y fl -x ll )a ll (y)+ d d y (1 - y^ + x^a^y), (3.1) 

where x now refers to the center of the link, i.e. X H IS cL half-integer. It is remarkable 
that this construction is gauge covariant. When one performs a gauge transformation 



= ~~ + $x-p,/2- (3.3) 



%(a;) = a^x) - d^x) (3.2) 
in the continuum, this induces a lattice gauge transformation 

The lattice gauge transformation is simply given by the integrated continuum trans- 
formation, i.e. 

$*= / d*y<p(y). (3.4) 

Jc x 

Note that it is sufficient to consider Abelian gauge transformations because we have 
not yet switched on interactions. A ^-function renormalization group transformation 
for the gauge field takes the form 

exp( — 5" [A]) = J T>a exp(— s[a]) 

x II s ( A * x - / ddy ( 1 + y»~ x »W(y) - / dd y ( 1 ~y^ + x »W(y)), 

(3.5) 

where s [a] is the continuum gauge action. In momentum space the lattice gauge 
field is given by 

A.(p) = V a.(p + 27r/)n,(p + 2tt/)(-1)^, IL.(p) = ^II(p). (3.6) 

lEZ d 

Because the link centers have a half-integer component x^, the gauge field is an- 
tiperiodic over the Brillouin zone in the /i-direction. For technical reasons we work 
in the Landau gauge = 0> i- e - the measure Da contains ^-functions 

3(^2/1 Pv a v(p))- Then the continuum action takes the form 

s[a] = J^yJ ddp ^(-p)^)" 1 ^)' = p" 2 ' ( 3 - 7 ) 

where S 9 (p) is the free gluon propagator in the Landau gauge. As in the fermionic 
case we want to optimize the renormalization group transformation such that the 
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fixed point action becomes as local as possible. Therefore we generalize the renor- 
malization group transformation of eq.(3.5) by smearing the ^-function. For this 
purpose we introduce an auxiliary lattice vector field D^, analogous to the aux- 
iliary fermion field rj. Also for we introduce a mass term ^aD^D^ to smear 
the ^-function. This alone would, however, not lead to an optimally local fixed 
point. Therefore we introduce in addition a kinetic term ^f> 2 D^D^ for the aux- 
iliary field. This will allow us to optimize the parameters in the renormalization 
group transformation such that the fixed point action becomes ultralocal even for 
2-d configurations. The introduction of a kinetic term for an auxiliary field was 
used before for staggered fermions [8, 25]. It is easy to show that the partition 
function remains invariant when one smears the ^-function in this way. We write 
the renormalization group transformation as 



exp(-S[A]) 



+ 



+ 



J VaVDexj>{- JL_ J d d p ^(-p)^)" 1 ^)} 
exp { - j^y d J d'p^D.i-p^a + jpDD.ip) 

iD^-p^A^p) - a »(P + 27r/)n^(p + 2tt/)(-1) 
iez d 

J VaVDexj>{- JL_ J d^-a^-p)^^ 1 a,{p) 
iD tl (-p)a tl (p)IL tl (p) } 

eXP {-(2^/ B 

>M-p)AM]}- 



} 



d d p 



(3.8) 



Performing the Gaussian integral over is equivalent to solving a classical equation 
of motion which yields 

<v(p) = iS^Tl^D.ip). (3.9) 
Since the continuum gauge field obeys the Landau gauge condition ?w(p) = 



we 



find 



Y,P.nv)^{p)D.{p) = Q. 



(3.10) 



Using p^IL^p) = p^IL^p) one obtains = 0> i- e - the auxiliary field 

obeys the standard lattice Landau gauge condition. Now we take the solution of 
eq.(3.9) and insert it in eq.(3.8). This is equivalent to integrating out the continuum 
gluons and yields 

exp(-S[A]) = JvD^ V [--^-y d Jj d p[iD.(-p)A.(p) 



+ ^(-P)E nv + 27r/)n,(p + 2tt/) 2 + a + IP^D.ip)] }. 



iez a 



(3.11) 
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Finally we integrate out the auxiliary vector field, i.e. we again solve a classical 
equation of motion, which now takes the form 

DM = -iA-Cp)- 1 ^), (3.12) 

with 

A£(p) = £ $ 9 (P + 27T0IW + 2tt/) 2 + a + 7 pJ. (3.13) 

iez d 

Combining eq.(3.12) with eq.(3.9) we obtain 

*M = ^mM&W-'Mp)- (3.14) 

This equation will be useful in sections 4 and 5 when we construct composite opera- 
tors and determine the vertex functions. The auxiliary vector field is in the standard 
lattice Landau gauge and hence one obtains 

5>A*(p)-%(p) = 0. (3.15) 

We call this the 'fixed point lattice Landau gauge' condition, because it is the one 
that arises from blocking gauge fields that obey the Landau gauge condition in 
the continuum. Then the measure DA contains ^-functions 6(^2 p^A^p) -1 A^p)). 
Inserting eq.(3.12) in eq.(3.11) one arrives at the fixed point action 

S[A] = j^y J B d d p l -A„{-p)Al{p)^A„{p), 

A£(p) = X(P + 27r/ )" 2 n^(p + 27r/) 2 + a + 7^- (3-16) 

iez d 

Here A^(p) is the lattice gluon propagator in the fixed point lattice Landau gauge. 
Of course, we would like to express the fixed point action in a gauge invariant 
way, i.e. we would like to go back to an unspecified gauge. Suppose the gauge 
field in an arbitrary gauge is given by ®Ap. Then ®Ap is related to the gauge 
field A^ in the fixed point lattice Landau gauge by a gauge transformation $, i.e. 
®Ap(p) = A^p) — ip^^p). The gauge condition (3.15) is an equation for $ which 
yields 

Inserting A^ into eq.(3.16) and dropping the index $ on ®Ap, one obtains the fixed 
point action in a gauge invariant form 

S[A] = j^y d J B d d p l -A„{-p)Al u {p)^A„{p), 
KM- 1 = A^)- 1 ^- ^ ^'f' (3-18) 
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Again we are looking for the optimally local fixed point action. It turns out 
that a and 7 can be optimized such that the action becomes an ultralocal single 
plaquette action for 2-d configurations. To see this it is advantageous to work with 
field strength variables F that live on the plaquettes. In the standard lattice Landau 
gauge the gauge potential in d = 2 is related to the plaquette variables by 

Mp) = ip- 2 e^P»F(p). (3.19) 

Inserting this into eq.(3.18), which is valid in any gauge — and hence also in the 
standard lattice Landau gauge we are using now — one finds after some manipula- 
tions 

S[F] = J^yJ B d2p I f (-p^ f ^ p^' 1 = p" a9 ^ + P2 Af (p). (3.20) 

Using eq.(3.13) one obtains 

P(P)~' = ^nP + 27r/)[P?n 2 (p + 27r/) 2 +^n 1 (p + 27r0 2 ] + a(p?+pD + 27P?P2 
;ez 2 

= ^n(p + 27r/) 4 + a(p 2 +p 2 ) + 27P 2 P2 



2 



(1 - lp 2 )(l -\p 2 2 ) + *(vl + f 2 ) + 2 7 plp 2 2 . (3.21) 



For 



a = -, 7 = , (3.22) 

6' ' 72' y J 

one arrives at p{jp) = 1. This corresponds to the standard plaquette action which 

is ultralocal. Hence at least in d = 2 the optimal choice for the smearing term is 

-^D^l — j^ffyD^. We use the same parameters in d = 4 and we write the fixed 

point action in coordinate space as 

S[A] = J] X -A^ xPllv {x - y)A„, y . (3.23) 

Performing a Fourier transformation of eq.(3.18) one finds the values of p^^z) that 
are given in table 3. As shown in fig. 3 even in d = 4 the fixed point action is 
very local (although not ultralocal). In [14] we will systematically vary a and 7 to 
optimize locality in d = 4. There a parametrization of the action in terms of closed 
loops will be given as well. 

Also for the gluons we would like to see that the spectrum of the lattice theory 
is perfect. We proceed in complete analogy to the quarks. For simplicity we work 
in the fixed point lattice Landau gauge. We construct an operator 

1 r 

Ai(p)x d = 7T dp d Ai(p)ex.p(ip d x d ), (3.24) 
2tt J-k 
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Z\ 


Z2 


Z 3 


z 4 


pll{ z ) 


Z\ 


Z2 


z 3 


z 4 
















4.1891 


0.5 


0.5 








-0.7120 





1 








-0.4437 


0.5 


0.5 


1 





-0.0640 





1 


1 





-0.1050 


0.5 


0.5 


1 


1 


-0.0106 





1 


1 


1 


-0.0324 


0.5 


0.5 


2 





0.0023 





2 








0.0030 


0.5 


0.5 


2 


1 


0.0001 





2 


1 





-0.0008 


1.5 


0.5 








0.0143 


1 











-0.0831 


1.5 


0.5 


1 





-0.0028 


1 


1 








0.0267 


1.5 


0.5 


1 


1 


-0.0009 


1 


1 


1 





-0.0038 


1.5 


1.5 








-0.0024 


1 


1 


1 


1 


-0.0033 


1.5 


1.5 


1 





0.0002 


1 


2 








-0.0024 


1.5 


1.5 


1 


1 


0.0002 


2 











0.0027 


1.5 


1.5 


2 





0.0001 


2 


1 








-0.0010 


2.5 


0.5 








-0.0005 


2 


1 


1 





0.0003 


2.5 


0.5 


1 





0.0001 


3 











-0.0001 


2.5 


1.5 








0.0001 



Table 3: Some values of pn(z) and p\2(z) in d = 4. The parameters a and 7 are 
chosen such that the action is ultralocal for 2-d configurations. Note that pn(z) 
is symmetric under permutations of z 2 , z 3 and Z4 while p\2(z) is symmetric under 
exchange of Z\ and z 2 as well as z 3 and Z4. 
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Z3 
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i1e-06 



1e-09 




Figure 3: Exponential decay of \p llll {i,i,i,i)\ for the free gluon (along the diagonal). 
The straight line is an exponential fit \p llll {i,i,i,i)\ oc exp(— ci) with c = 6.937 
excluding the point (0,0,0,0). 

and we consider the 2-point correlation function 

i r 

{Ai,(2,o)Ai(p)x d ) = 7T I dpd Af(p)exp(ip t/ a; t/ ) 

1 Z" 7 

— / dp d [ V(p + 27r/)- 2 n,(p + 2tt/) 2 + a + 7 p 2 ] exp(ip d x d ) 

1 , , ^ n 



iez d 



= j_ r d v — 

ZW-oo (p + 27rr) 2 +p 2 Vp 8 + 27r/ 8 7 11 V Pj + 2^ 

(Pc/ \ 2 
— exp^p^rf) + (a + J$)8 Xd ,o 

= J] C'(p + 27rr)exp(-E(p + 27r/> ti ) + ( a + 7 p 2 )^,o. (3.25) 
fez**- 1 

Again the energies of the lattice theory are given by the locations of the poles. Now 

E(p + 27rf) 2 = -p 2 = (p + 27rf) 2 , (3.26) 

which is the energy of a massless gluon with momentum p + 27r/, i.e. again the 
spectrum is identical with the one of the continuum theory. 

As for the quarks we also want to discuss renormalization group transformations 
with a finite blocking factor n. In this transformation we construct blocked link 
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variables 



— V- Va 



A' 



c-\-fi/2-\-ifi ' 



(3.27) 



We average over all straight line connections between corresponding sites in adjacent 
blocks. The factor j3 n is a wave function renormalization. In momentum space the 
renormalization group transformation takes the form 

eMS'lA']) = J VAVD ^v{~j^yJ B d d P ^(-p^p)" 1 ^)} 
X 6XP { " ( 2 7) d ^,^[^ ( " p)(a " +7 "^ ) ^ (p) 

+ iD,(-p)lK(p) " & £ MP + — + — )(-i)'J-]l }• 



(3.28) 



The smearing parameters are now given by a n and 7„. Also n n/J (p) = 11^ (rap) / 11^ (p) 
and the summation extends over integer vectors with components l'^ £ {1, 2, 
Integrating out the fine lattice fields one obtains the coarse gluon propagator 



TV 



n 



n 



+ a + 4 7 sin 2 (^^)]n„,(^±^) 2 + a „ + 7 „^ 



n 



+ 



d-2 



n 



£ (p, + 27T/;)- 2 n,(p + 27T/') 2 



6ra 2 



A 2 \ i A 2 i , A 2 

Pn) + 1rfPn + (X n+lnP ll - 



(3.29) 



This is the perfect propagator for blocking from the continuum provided that j3 n = 

n (d-2)/2^ a j n 2 _|_ a ^ _ a an( J -y/ n 4 _|_ ^ _ 0,^2 _ ^^5^4 _ ^_ fp^g va } ue Q f ^ 

consistent with dimensional considerations because the dimension of a gauge field is 
(d — 2)/2. One also obtains 



OLr 



a 1 



7n 



7 1 



7r 



6ra 2 



7JT 



(3.30) 
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In particular = a and 7^ = 7 because blocking from the continuum corresponds 
to an infinite blocking factor. Using the optimized values a = 1/6 and 7 = —1/72 
obtained in eq.(3.22) for blocking from the continuum one finds for general block 
factor n 




(3.31) 



In particular for block factor 2 the values are a 2 = 1/8 and 72 = —1/128. 



4 Classically perfect fields, composite operators 
and currents 

In this section we construct classically perfect expressions for quark and gluon fields, 
which we use to construct lattice versions of quark-gluon composite operators of the 
continuum theory. These operators may serve as a first approximation to quantum 
perfect operators that one would like to use e.g. in numerical simulations of hadronic 
matrix elements. In particular we construct the Polyakov loop as well as current 
operators. In general, the correlation functions of the classically perfect fields are 
not completely free of cut-off effects. However, as we show explicitly for the quark 
2-point function and for the correlation function of two Polyakov loops, the cut-off 
effects are exponentially suppressed. We also construct other lattice currents which 
have an exact interpretation in terms of continuum physics. They are useful when 
one wants to discuss Ward identities or investigate anomalies in the perfect action 
approach. For example, in [22] it has been shown that the exact lattice axial current 
has the correct axial anomaly in the 2-d lattice Schwinger model, whereas the naive 
Noether current is conserved. 

Let us consider the continuum quark field ij). As it was observed before in scalar 
field theory [10], at the classical level the continuum field can be reconstructed from 
the lattice field This was achieved in eq.(2.11) were we found 

M-p) = ^(-pJA^p)-^^)^^), MP) = 6 f (p)IL(p)Af(p)- 1 ^(p). (4.1) 

The correlation function (ij) c {x)il) c {y')) is invariant only against simultaneous trans- 
lations of x and y by lattice vectors, although x and y may be arbitrary points 
in the continuum. To restore full continuum translation invariance we consider 
J d d z (ip c (x + z)ip c (y + z)), i.e. we average the correlation function over a hyper- 
cube c centered at the origin. Note that the displacement between the quark fields 
is fixed. In momentum space this corresponds to 

I(&(-pty c (p)) = ^fn^A^)- 1 . (4.2) 

Here V is the space-time volume. Obviously the result differs from the correct 
continuum expression y(ip(— p)ip(p)) = S^(p). The same is true in scalar field theory. 
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Therefore the 2-point function of the classically perfect field is not perfect — it is 
different from the one of the continuum theory, and hence it is not free of cut-off 
effects. Still the classically perfect field is an interesting object, because — as we 
will see now — the cut-off effects in its 2-point function are exponentially suppressed 
at large distance. From eq.(4.2) we infer 

(M-P)MP)) = m-p)Tp(p)) S^pMpYA^p)- 1 . (4.3) 

In coordinate space this implies that J d d z(ip c (x + z)ip c (y + z)) is a convolu- 
tion of the continuum 2-point function {^{x^ijj^y)) with the Fourier transform of 
6? (p)Il(p) 2 A-^(p) -1 which is defined at any point in the continuum. We want to 
show that this function decays exponentially. For this purpose we restrict it to the 
lattice points. In momentum space this corresponds to p £ B and a summation over 
integer vectors 

^^ / (p+27r/)n(p+27r/) 2 A / (p+27r/)- 1 = [A f (p)-a]A f (p)- 1 = 1-oA^p)" 1 . (4.4) 
iez d 

Here we have used eq.(2.10) as well as the periodicity of A-^(p). In coordinate space 
the result corresponds to S XtV — a [7 ^{x — y) + A(a; — y)]. As we have seen in section 
2 this function does indeed decay exponentially. Hence the 2-point function of the 
classically perfect field differs from the correct 2-point function of the continuum 
theory only by exponentially small cut-off effects. The same was concluded in the 
first paper of ref.[10] for scalar field theory. 

Of course we can also construct the classically perfect gluon field. In fact, this 
has already been done in eq.(3.14) where we found 

<V = S^p^p^pY'A^p). (4.5) 

Note that in this equation is in the Landau gauge and is in the fixed point 
lattice Landau gauge. Using eqs.(4.1) and (4.5) one can construct a classical lattice 
version of any quark-gluon composite operator 0[^,-^,a] of the continuum theory 
simply by identifying 

c [tf,*,A] = 0[i,i,flj. (4.6) 

The correlation functions of the composite operators will not be perfect. However, 
as for the 2-point function of the classically perfect quark field one expects only 
exponentially small cut-off effects. This may be useful e.g. in numerical simulations 
of hadronic matrix elements. As an example of a composite operator we consider 
the Polyakov loop 

4>(x) = J dx d a d (x,x d ). (4.7) 

Here we restrict ourselves to the leading order in the gauge coupling. Then path 
ordering problems do not occur. In momentum space we have (j)(p) = a d (p,0). The 
Polyakov loop expressed in terms of lattice variables is then given by 

Mp) = VMp, o)A*(p, o)- 1 ^(p, 0). (4.8) 
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In general the classically perfect Polyakov loop is defined at any point in the con- 
tinuum coordinate space. If we restrict ourselves to lattice sites x we can write 

= <t>c{x) = ^u(x- y)$g, Qjj = ^ A dUyd) , (4.9) 
y Vd 

where $^ with y £ Z rf_1 is the standard lattice Polyakov loop. In momentum space 
the above kernel is given by 

"(P) = Yl + ^)~ 2 MP + 2ttI,0)A%p,0)-\ (4.10) 



iez d 



where 



d-i 



d-i , 2 



£<?+M- 2 n(diri) +«■ < 4 - n > 



,- €Z d-i ,=i ^ 



First we consider d = 2 and find 

P? , .._ 1 1 

f?5i - 

/i£Z 



^°) = i;5rTU? +a = ^-5 +a ' (4 - 12) 



i.e. for the value a = 1/6, which yields an ultralocal fixed point action in d = 2, 
one obtains A^pijO) -1 = j>\. Inserting this in eq.(4.10) one finds u;(0) = 3/4 and 
a>(±l) = 1/8. At all larger distances uj(z) vanishes. Hence the classically perfect 
Polyakov loop in d = 2 is ultralocal and is given by 

$ c ,x = ^($x+i + 6$. + $.-1) (4.13) 
o 

Some values of uj(z) for d = 4 are given in table 4 and in fig. 4. Again the classically 
perfect Polyakov loop is very local (but no longer ultralocal). 

Next we consider the correlation function of two Polyakov loops. It takes the 
form 

^(M-P)MP)) = ^Ofn^OfA^O)- 1 . (4.14) 

Again this is different from the continuum result y (</)(— p)</>(p)) = S 9 (p,0). As 
for the quark 2-point function we like to show that the cut-off effects are expo- 
nentially small. In coordinate space the left hand side of eq.(4.14) corresponds to 
Jf o d d ~ 1 z((/) c (x + z)(/) c (y + z)), where / is a (d— l)-dimensional cube centered at the 
origin. The correlation function averaged over / is invariant against arbitrary con- 
tinuum translations. It is a convolution of the correct continuum result ((j)(x)(j)(y)) 
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^2 


^3 


cj(z) 











0.48709 


1 








0.05188 


1 


1 





0.01354 


1 


1 


1 


0.00460 


2 








0.00015 


2 


1 





0.00010 


2 


1 


1 


0.00004 


2 


2 





-0.00007 


2 


2 


1 


-0.00004 


2 


2 


2 


-0.00002 


3 








0.00001 



Table 4: The largest values ofuj(z). We have chosen a = | such that the classically 
perfect Polyakov loop is ultralocal in d = 2. Note the symmetry under permutations 
of the Zi. 

with the Fourier transform of S g (p, 0)I[d(p, 0) 2 A^(p, 0) _1 . Again we restrict this 
function to the lattice points which in momentum space leads to 

J] S 9 (p + 2tvI, 0)U d (p + 2tvI, 0) 2 A%p + 2tvI, 0)" 1 = 

fez**- 1 

[A%p,0) - aJASfoO)- 1 = 1 - aASfoO)" 1 . (4.15) 

Because p^(x — y) decays exponentially (as we saw in section 3) the resulting cut-off 
effects are again exponentially small. This is in agreement with a result obtained 
in [10] for a different renormalization group transformation. From eq.(4.14) one 
obtains an expression for the static quark- antiquark potential 

V{f) = -J^yZT J dd ~ X V 6 9 (P, 0) 2 n d (p, 0) 2 A%p, 0)" 1 exp(ijf ■ f). (4.16) 

Note that V(—r) = V(r). First we consider the potential in d = 2. There the 
integral can be performed analytically and the potential takes the form 

( r/2 for r > 2 

V(r) = I -r 5 /120 + r 4 /12-r 3 /3 + 2r 2 /3-r/6 + 4/15 for l<r<2 (4.17) 
{ r 5 /40 -r 4 /12 + r 2 /3 + 7/30 for < r < 1, 

which is four times continuously differentiable. We have fixed an arbitrary constant 
in the potential such that it agrees with the continuum result r/2 at large distances. 
It is remarkable that the potential obtained from the classically perfect Polyakov 
loop then agrees exactly with the continuum potential for r > 2. The potential 
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1 2.3 4 

i 

Figure 4: Exponential decay of \cj(i,i,i)\ for the classically perfect Polyakov loop 
(along the diagonal). The straight line is an exponential fit \cj(i,i,i)\ oc exp(— ci) 
with c = 5.646 excluding the point (0,0,0). 

in d = 2 is shown in fig. 5. Even for 1 < r < 2 — and hence at any distance 
larger than one lattice spacing — it differs from the continuum result by less than 2 
percent. It should be noted that the potential of the standard Polyakov loop agrees 
with the continuum result. However, it is restricted to lattice points, whereas the 
potential of the classically perfect Polyakov loop can be evaluated at any distance. 
The potential in d = 4 is compared to the continuum potential — l/47rr in fig. 6. We 
show the potential derived from the classically perfect Polyakov loop both along the 
space diagonal and along an axis. The two agree at the 1 percent level at any distance 
r, i.e. rotation invariance is restored with good accuracy. This is in contrast to the 
standard Wilson potential. Even at r = 3 its values at f = (3,0,0) and (2,2,1) 
differ by more than 5 percent. The potential of the classically perfect Polyakov 
loop along the diagonal, i.e. V(r / \/3,r / \/3,r / deviates from the continuum 
potential by less than 0.4 percent at any distance larger than 1.5 lattice spacings. 
In [10] the potential of the classically perfect Polyakov loop has been determined 
for a different renormalization group transformation at distances that correspond 
to lattice vectors. Even at r = (1,0,0) it agrees with the continuum answer at the 
percent level. This is not the case for our renormalization group transformation 
when we put a = 1/6 providing ultralocality in d = 2. However, as an alternative 
criterion one may choose a such that V(r) is as close as possible to the continuum 
potential. This would suggest a 1/13. For that value the agreement with the 
continuum potential is of the same quality as in [10]. In [14] we will systematically 
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Figure 5: The static quark- antiquark potential in d = 2 obtained from the classically 
perfect Polyakov loop (solid curve) compared to the continuum potential (dashed 
line). 

vary a and 7 to optimize both the potential and the locality of the action in d = 4. 

Next we turn to the currents. We start in the continuum and consider the current 

j ll (x)=^(x)T^(x), (4.18) 

where Y ^ = 7^ for the vector current, and Y ^ = 7^75 for the axial vector current. 
Of course, more general expressions for Y ^ are also possible; for example Y ^ can be 
a flavor-nonsinglet current as well. In momentum space eq.(4.18) takes the form 

We define lattice currents by inserting the classically perfect fields of eq.(4.1) into 
eq.(4.19) and we obtain 

^c(p) = tAz / d d q^(p- q )Af(p-q)- 1 V 6f(p-q-2Tvm)Y„ 

raeZ" 

X S f (q + 2Tvm)YL(p- q-2Tvm)YL(q + 2Tvm)A f (q)V(q). (4.20) 

We also want to define currents that are genuine lattice objects. Still they should 
have an exact interpretation in terms of continuum physics. A natural definition is 
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Figure 6: The static quark- antiquark potential in d = 4. The dashed curve is the 
continuum Coulomb potential —l/Anr. The potential obtained from the classically 
perfect Polyakov loop along the space diagonal and along an axis is represented by the 
solid and dotted curve, respectively. The potential can be evaluated at any distance, 
whereas the potential of the standard Wilson theory, represented by the diamonds, is 
limited to lattice vectors. 

the total flux of the continuum current through a (d — l)-dimensional face jX that 
separates two hypercubes c x _f 1 / 2 and c x j r f 1 /2. The corresponding lattice current 
associated with the link connecting the lattice points x — fi/2 and x + fi/2 is given 
by 

U = f d^yj^y). (4.21) 

We call this the 'exact lattice current'. Its standard lattice divergence is by Gauss' 
law identical with the divergence of the continuum current integrated over a hyper- 
cube c x 

SJ X = ^(<W,V2 - Jn,x-fi/2) = / d d y d^j^y). (4.22) 
In momentum space the exact lattice current takes the form 

Up) = Y, Up + 27r/)n^(p + 2tt/)(-i)^, il^( p ) = ^n( P ). (4.23) 

H P» 

The current is antiperiodic over the Brillouin zone in the /i-direction because it lives 
on the link centers, which have a half-integer component x^. As it stands the lattice 
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current is expressed in terms of continuum fields. We now express it in terms of 
the lattice quark fields by using eq.(4.20) 

X ^ S f (p + 2Tvl-q-2Tvm)T tJ S f (q + 2Tvm)A f (q)- 1 ^(q) 

l,mEZ d 

)n^(p + 27r/)(-l)^. (4.24) 

The exact lattice current is analogous to the lattice field $ because it is a genuine 
lattice object that still has an exact interpretation in terms of continuum physics. In 
particular its lattice divergence is the divergence of the continuum current averaged 
over a hypercube. This property was essential in a study of the axial anomaly in 
the Schwinger model [22]. It should be emphasized that correlation functions of the 
exact currents are not free of cut-off effects. In fact, their cut-off effects are not even 
exponentially small. Still the exact current is a natural object with useful properties 
in the perfect action approach to lattice field theory. 

5 The perfect quark-gluon and 3-gluon vertex 

In this section we switch on the interactions between quarks and gluons. However, 
we limit ourselves to the leading order in the gauge coupling e. Now the non- 
Abelian nature of the fields becomes important and we no longer suppress the color 
indices. The continuum quark field is denoted by ijj 1 with i £ {1,2,..., TV} and 
the continuum SU(N) gauge field is denoted by with a £ {1,2, ...,7V 2 — 1}. The 
gluon field strength is given by 

= - d v al{x) + ier bc al(x)al(x). (5.1) 

Here f abc are the structure constants of the SU(N) algebra spanned by the Her- 
mitean generators A a , such that [A a ,A b ] = 2if abc X c and Tr(A a A b ) = 2S ab . The con- 
tinuum QCD action is given by 

s[M,a] = J d d x{^(x)f^(x) + nx)[^(d^ (5.2) 

In momentum space the resulting quark-gluon vertex takes the form 

v[M,a] = j^- d J d d pd d q ^(-p)i^a;(p - q)X^(q), (5.3) 

and the 3-gluon vertex is given by 

v[a] = jdr J ddpddq ^ abc <(-*" - - «K(«)- ( 5 - 4 ) 
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Again we want to integrate out the continuum fields and derive the perfect quark- 
gluon and 3-gluon vertex of the lattice theory. Before we can do so we need to define 
the renormalization group transformation up to 0(e). This can be done in many 
ways, and we do not want to make a specific choice yet. In general, for the quarks 
one can write 

= ^ i>\p + 27r/)n(p + 2tt/) 

iez d 

+ 6 X) T^y I d<lq K ^ P + 2?r/ ' q + 2 * l ) a > ~ 9)^-^(9 + 2tt/). (5.5) 

Here K li {jp,q) is a regular kernel that specifies the renormalization group transfor- 
mation. It describes how the continuum quark field in a hypercube c x is transported 
to the lattice point x. For a certain class of renormalization group transformations 
the kernel obeys an additional gauge covariance constraint. In these cases one re- 
quires that a continuum gauge transformation performed on ij) and induces a 
lattice gauge transformation of This implies the constraint (p^ — q^K^p^q) = 
II(p — q)IL(q) — II(p). It is not really necessary to impose this constraint. In fact, in 
the next section we will select a specific renormalization group transformation for 
QCD for which the kernel is identically zero. To determine the kernel that corre- 
sponds to a given renormalization group transformation with a finite blocking factor 
n one proceeds as follows. First one blocks from the continuum to a lattice of spac- 
ing 1 and performs one additional step of the block factor n transformation. Then 
one derives a recursion relation for the kernel by demanding that the result agrees 
with the one obtained by blocking from the continuum to a lattice of spacing n. 

Compared to the free fermion calculation of section 2 we now have an extra term 

T \MM = J^ d Jd d pd d qm-p)K„(p,q)a;(p-q)\^(q) 

- ^(-p)K ll (-q,-p)al(p-q)X^(q)] (5.6) 

in the renormalization group transformation and hence in the exponential of eq.(2.5). 
The additional terms v[ij),il),a\ and T[^,-^,a] modify the exponent E[i/>,ij>] + E[a] 
of the free theory to E[i/>,ij>, a] = E[ij),il)] + E[a] + v^,^, a] + T[^,-^, a]. Since the 
additional terms are 0(e) they change the equations of motion (and hence their 
solutions) to ij) c + 6ij> c and ij) c + 6ij> c and a c + Sa c , where 6ij> c , 6ij> c and Sa c are of 0(e). 
Inserting the new solution into the new exponent yields 

E[$ c + 6ij> c , i> c + 6ij> c , a c + Sa c ] = 

E[^ c + ^ c ,i + 6i}) c } + E[a c + Sa c ] + v[4> c ,tl> c ,a c ] + T[^ c , ij) c , a c ] + 0(e 2 ) = 

+ E[a c ] + v[ip c ,ip c ,a c ] + T[ip c ,ip c ,a c ] + 0(e 2 ) = 
S[V, *] + S[A] + V[V, *, A] + 0(e 2 ). (5.7) 

The essential observation is that ij) c , ij) c and a c are solutions of classical equations of 
motion, i.e. extrema of the corresponding exponents. This implies E[ij> c + ^ c ,-^ c + 
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6i/> c ] = Eli/jci/jc] + 0(e 2 ) and E[a c + Sa c ] = E[a c ] + 0(e 2 ). Hence it is sufficient to 
insert the old solutions of the equations of motion into the additional terms. This 
allows us to read off the perfect quark-gluon vertex 

V[V,V,A] = v[i> c ,i> c ,a c ] + T[i> c ,i> c ,a c ] 

= j^r L ddpddq ** ( " p)e ^ (p ' q)A » {p ~ «) ( 5 - 8 ) 

Using eqs.(2.9), (2.11) and (3.14) one identifies the quark-gluon vertex function as 
V.(p,q) = A^pr'AKp-q)- 1 Y, S°(p + 2Tvl-q-2Tvm) 

l,mEZ d 

X n^(p + 27r/ - ? - 27rm)(-l/" +m " 

X S f (p + 2Trl)i'-f ll S f (q + 27rm)II(p + 27r/)II(g + 2?rm) 

+ K^p + 2?r/, q + 2-wm)6 i (q + 27rm)II(g + 2?rm) 

- K tl (-q-2Tvm,-p-2Tvl)S f (p + 2Tvl)U(p + 2Tvl) A f (q)-\ (5.9) 

This expression is in the fixed point lattice Landau gauge. To transform it to a 
general gauge one proceeds in complete analogy to the free gluon case. We will do 
this explicitly in the next section when we choose a specific renormalization group 
transformation. 

Also the renormalization group transformation for the gluons needs to be defined 
up to 0(e). We write 

KiP) = E* + 27r/ )n,(p + 27rZ)(-iy<< 
iez d 

+ e E TTw / ^ f abcR ^ P (P + 27r/ > 9 + 27r 
^ (2?r) J 

X ai(p-q)a c p (q + 2irl)(-l) 1 ". (5.10) 
Compared to the free gluon case of section 3 we now have an extra term 

T[a] = (2^ / ddpddq f abCD :(-p) K ^ M(p - <?K(<?) (5-n) 

in the exponential of the renormalization group transformation eq.(3.8). In complete 
analogy to the quark-gluon vertex one obtains for the 3-gluon vertex 

V[A] = v[a c ] + T[a c ] 

= T^XmI d d pd d qr bc Al(-p)eV t ^ p (p,q)A b l/ (p-q)Al(q). (5.12) 
l Z7r J Jb 2 
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Using eq.(3.14) one identifies 

V„„ p (p, q) = A^p)- 1 Al(p - qV^qV £ 6°(p + 2tt/ - q - 27rm) 

l,mEZ d 

X n„(p + 2tt/ - q - 2Trm)(-l) lv+mv 

X £ 9 (p + 2?r/)[p^ + 2tt/^ - q^ - 2-Km ll ]8 up 8 9 (q + 2?rm) 

x n^(p + 27r/)n p ( ? + 27rm)(-l)^ +ra - 

+ K^ p (p + 2tt/, ? + 27rm)^( ? + 2Tvm)U p (q + 27rm)(-l) ra '] .(5.13) 
Also this expression is in the fixed point lattice Landau gauge. 



6 A nonperturbative version of the renormaliza- 
tion group transformation 

Before one uses perfect actions in numerical simulations of light quarks one would 
also like to incorporate nonperturbative effects. This requires numerical work along 
the lines of [4, 10]. Then it is essential to use a renormalization group transformation 
with a finite blocking factor n. In practice one would use n = 2. Below we work with 
general n. Our goal is to formulate a nonperturbative version of the renormalization 
group transformation whose perfect action reduces to the one of blocking from the 
continuum at least for weak fields. Then we know that the perfect action is very 
local (in lower dimensions even ultralocal), and therefore well suited for numerical 
simulations. 

The nonperturbative formulation of lattice gauge theories works with compact 
link variables U^^ (E SU(N) connecting neighboring lattice points x and x + fi. Note 
that here we have not associated x with the link center; x^ is now an integer. The 
parallel transporters are related to the gauge potential by 

= exp(^ + - /2 A a ). (6.1) 

As usual we have rescaled eA^ — > A^. Under a gauge transformation g the link 
variables transform as 

9 U^ x =g x U^ x g + x+[l . (6.2) 

In the renormalization group transformation we want to couple a blocked link U' x , 
on the coarse lattice to an average of n d parallel transporters on the fine lattice. The 
blocked link connects the centers of two neighboring n d blocks on the fine lattice. 
For each point x £ x' in the first block we construct a parallel transporter as a 
product of n consecutive links starting at x and ending at the corresponding lattice 
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point x + n\L in the second block. Then we couple the parallel transporters to the 
blocked link U' x , by constructing a traceless anti-Hermitean matrix 

^ n — 1 

B »y = - d J2 lo ^KAU (6.3) 

for each blocked link. A gauge transformation g' on the coarse lattice 

'Ks=&Ks&+nfi + (6.4) 

induces a gauge transformation g x = g' x , for x £ x' on the fine lattice. Hence the 
anti-Hermitean matrix transforms as 
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B^y = g x ,B^ x ,g x , . (6.5) 

Note that transforms with the same g' on both sides. Hence it is a charged vector 
field in the adjoint representation, not a lattice gauge field. In the renormalization 
group transformation we will need the covariant second derivative of in the fi- 
direction 

A^U'jB^y = 1B^ X , - U'^ x ,B^ x i +nii U'^ xl + - U^_ nfi + B^ x i- nii U'^ xl _ nii . (6.6) 

(There is no summation over fi.) Note that the fields are shifted by n(l because the 
blocked lattice has lattice spacing n. Under a gauge transformation g' the covariant 
second derivative transforms as 

A^'ur'B.y = g^A^B.yg'^. (6.7) 

Now we are ready to formulate the nonperturbative version of the renormalization 
group transformation for QCD. The effective action for the coarse variables after a 
factor n renormalization group step is given by 

exp(-5'[*', *',[/']) = J exp(-S[*, *,[/]) f[U] 

x' x£x' x£x' 

x exp{ -^Y / ^[KA^+7nA.[U'])- 1 B^ xl ]}.(6.S) 

X 1 ,[1 

The functional f[U] ensures the proper normalization of the partition function and 
is given by 

fill]- 1 = J W'exp{ --L^Tr^K + T^tC/'])" 1 ^]}. (6.9) 
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For weak gauge fields the above transformation reduces to the discrete transforma- 
tions of sections 2 and 3. Hence we know that the parameters 




(6.10) 



give a very local perfect action. From eq.(6.8) one automatically obtains j3 n = n 2 . 
(/3„ was introduced in eq.(3.27)). A fixed point is reached only for f3 n = n^ d ~ 2 ^ 2 
which limits us to d = 4. This is no surprise because only in d = 4 non-Abelian 
gauge fields are asymptotically free. Then the fixed point of the interacting theory 
is close to the Gaussian fixed point of the free theory. 

It may seem that the above transformation is not gauge covariant, e.g. because 
the fine lattice quark field ty x is not properly parallel transported to the block center 
x'. This is, however, not required for the gauge invariance of the blocked action. 
In fact, performing a gauge transformation g' on the coarse lattice and the induced 
gauge transformation g with g x = g' x , for x £ x' on the fine lattice leads to 

x£x' x£x' 

(*^ + -^£ **tiM>K> ~ ~ d £ 9***) = 

n a z — ' n a z — ' 

x£x' x£x' 

» - ^ £ *-)(^' - ^ £ *-)■ (e- 11 ) 

n a A — 4 n a A — 4 

x£x' x£x' 

Using gauge invariance of the fine action and of the measure one observes that the 
coarse action is indeed gauge invariant. In particular, for the above transformation 
the kernel K^p, q) of eq.(5.5) vanishes identically and the vertex function of eq.(5.9) 
reduces to 

V„(p,q) = A'Cp^A-Cp - q)- 1 £ S°(p + 2irl-q-2irm) 

l,meZ d 

X n^(p + 27r/- q-2Tvm)(-l) l » +m » 

x S f (p + 2*1)^^ (q + 27rm)n(p + 2Tvl)U(q + 27rm)A / ( ? )- 1 .(6.12) 

Here we are in the fixed point lattice Landau gauge. To go back to a general gauge 
we perform the gauge transformation of eq.(3.17) both on the fermions and on the 
gauge fields and obtain 

V.(p,q) = ^(pr'A^p-q)- 1 £ S 3 (p + 2irl-q-2irm) 

l,mEZ d 

X H v (p + 2irl - q - 2Trm)(-l) lv+mv 
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X 



S f (p + 2Tvl)ijJ f (q + 27rm)n(p + 2Tvl)U(q + 2Tvm)A f (q)' 1 



E p (p - <i) 2 P K(p - q)' 1 

This expression contains products of up to five 7-matrices. They may be decomposed 
into 1, 7^, cr^, 75 and 7^75. The cr^-term represents a perfect 'clover' term. This 
decomposition will be carried out in [14]. There we will also parametrize the vertex 
function in terms of parallel transporters connecting the quark fields. Here we 
consider the vertex function for 2-d Abelian gauge fields. In addition we assume that 
the fermion fields are constant in the 2-direction, i.e. ^( Xl ,x 2 ) = ^(^1,0)5 ^( Xl , X2 ) = 
$( Ib o). Hence they are effectively 1-dimensional. Then the vertex function can be 
computed analytically in coordinate space, and to 0(e) the total action takes the 
form 

s[*,*,a] = Y,\f 2 x 

X 

X 

+e ra i^ " + ieA ltX+i/2 )* x+i - Ut x+i {l - ieA ltX+i/2 )* x } 

+ ("?) XI {**°2X(:ci - 2/1,2/1 - z 1 )ieA 2 ,yV z + ^ x a 3 6(x 1 - y 1 ,y 1 - zi)ei^$ z |. 

(6.14) 

Here 

F x = A l x _2/2 + ^-2,a;+i/2 ~~ A l,x+2/2 ~~ ^2,ai-i/2 (6.15) 

is the field strength living at the plaquette centers. Note that in two dimensions 

i 

v^F^ = v^F^ = -a 3 F. (6.16) 

Hence the last term in eq.(6.14) represents a perfect clover term. The nonvanishing 
values of x( x iiVi) an( i ^( x 1tUi) are given in tables 5 and 6, and their relative weights 
are illustrated in figs. 7 and 8. The vertex functions are peaked at X\ = —y\, i.e. 
where fermion and antifermion are at the same position. For fixed separation X\ +2/1 
the maximum occurs when the gauge field is close to the fermions, i.e. when X\ — y\ is 
small. It is remarkable that the vertex function is again ultralocal for 1-d fermions 
coupled to 2-d gauge fields, i.e. the coupling vanishes identically beyond short 
distances. The full 4-d vertex function will of course not be ultralocal. However, 
based on the experience with the propagators we expect it to be still very local. 

By construction, in the classical limit e — > the full renormalization group 
transformation reduces to the quadratic problem of section 3 even for arbitrarily 
strong abelian gauge fields. Since we know that for 2-d configurations the fixed 



32 



Xl 


2/i 




-2 


1 


1/80 


-2 


2 


1/120 


-1 





17/240 


-1 


1 


1/6 








19/60 



Table 5: The nonvanishing values of the effectively 1- dimensional vertex function 
Xixr.yr). Note that x(xi,yi) = x(~Xi,-yi) = x(yi,Xi). 





2/i 


0( x u 2/i) 


-1.5 


0.5 


1/80 


-1.5 


1.5 


1/120 


-0.5 


-0.5 


7/120 


-0.5 


0.5 


19/120 



Table 6: The nonvanishing values of the effectively 1- dimensional vertex function 
8(xi,yi) entering the perfect clover term. Note that 8(xi,yi) = 8(—Xi,—yi) = 
<%i, ap- 
point action is ultralocal at the quadratic level, it follows that for any 2-d abelian 
gauge field configuration the single plaquette Manton action [15] is classically perfect. 
This is a remarkable result, because it is valid for strong fields. Our construction of 
the renormalization group transformation appears very natural. The perfect action 
of the transformation reduces to the one of blocking from the continuum for weak 
fields. Its parameters can be chosen such that the perfect action is very local. Besides 
these nice analytic properties the transformation is well suited for numerical work. 
For example, in the classical limit e — > (and for the values of a n and 7„ given 
in eq.(3.31)) it turns out that f[U] = 1. This is very useful in the nonperturbative 
numerical determination of the classically perfect action beyond the weak field limit. 
This will be discussed in a subsequent publication [14]. 

7 Conclusions 

In this paper we have performed the analytic calculations that are necessary to 
construct a perfect action as well as perfect classical fields and their composite 
operators for full QCD. The technique of blocking from the continuum allowed us 
to derive closed analytic expressions for the perfect action up to O(e). This is not 
the case for other renormalization group transformations. In general, one can only 
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2.5 

1.5 

0.5 
-0.5 
-1.5 
-2.5 

-2.5 -1.5 -0.5 0.5 1.5 2.5 

Figure 7: Illustration of the relative weights of the nonvanishing values of the ef- 
fectively 1- dimensional vertex function x{_ x i,yi)- Note the peak at X\ = —y\ where 
fermion and antifermion are at the same position. 

derive recursion relations for the perfect propagators and vertex functions, but one 
cannot iterate to the fixed point analytically. We are convinced that it is very 
useful to have closed analytic expressions for the fixed point quantities, not only 
because they provide physical insight, but also because they simplify the numerical 
calculations that are necessary to turn the analytic results into a practical lattice 
action that can be simulated. 

An alternative approach could start from the fixed point of free staggered fermi- 
ons which has been constructed in [8, 25]. There one is limited to odd blocking 
factors, and one would use n = 3 [26]. Then also the blocking of the gauge field 
must be modified in a suitable way. Again one can work out the analytic expressions 
for the perturbatively perfect action. This is presently in progress [27]. 

For simulations it is essential that the renormalization group transformation 
we have used has a very local perfect action, provided its parameters are tuned 
appropriately. In particular, in lower dimensions (1-d for the fermions, 2-d for the 
gauge field) the perfect action of the free theory reduces to the ultralocal standard 
action. It is remarkable that even for arbitrarily strong 2-d abelian gauge fields 
the single plaquette Manton action is classically perfect. We expect that the fully 
nonperturbative perfect action of our renormalization group transformation will also 
be very local. Hence it should be possible to approximate it by a practically perfect 
action that can be simulated numerically. 

The obvious next step is to extend our results nonperturbatively. This requires 
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-2 
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Figure 8: Illustration of the relative weights of the nonvanishing values of the effec- 
tively 1- dimensional vertex function 8(xi,yi). Note that X\ + yi is the separation 
between fermion and antifermion. 

numerical work along the lines of [4, 10] which is presently in progress [14]. It is cer- 
tainly nontrivial to include quarks in the nonperturbative perfect action. A natural 
first step will be to consider heavy quarks. First, their cut-off effects are large and 
therefore a perfect action should be much better than the standard Wilson action. 
Second, they are easier to simulate than light quarks, and third, the gauge inter- 
action is weak at the scale of the heavy fermion mass such that our perturbatively 
perfect action may already be useful. In addition the perfect action of a heavy quark 
is more local than the one of a light quark. The ultimate goal is, of course, to include 
light fermions as well. 
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